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1 Introduction 

The communication between different brain struc¬ 
tures is crucial in cognitive functioning. In human 
subjects, such interactions have been studied nonin- 
vasively only by analyizing raw MEG and EEG sen¬ 
sor signals [1]. The conclusions resulting from these 
analyses are hampered by the fact that the sensitivity 
profiles of the sensors largely overlap. Furthermore, 
the emphasis has been on establishing relations be¬ 
tween the activities picked up by the sensors, and no 
attempts have been made to infer “causal” relations. 
The temporal correlation pattern between signals and 
phase differences might provide clues as to which of 
the activities leads with respect to the others, which in 
turn might provide a basis for the inference of cause 
and effect in brain activity. 

Electromagnetic source analysis overcomes the prob¬ 
lem of the overlapping sensor sensitivity profiles, 
making it possible to obtain more conclusive evidence 
for the origins of the observed activity. 

It therefore seems desirable to generalize the spa- 
tiotemporal source analysis procedure to incorporate 
structural relations between the activities of different 
sources. 

Interactions between sources may be conceptualized 
by the transfer of activity from one source to another. 
The interactions between cortical areas may then be 
analyzed by estimating the transfer function from one 
source to another. 

A major problem in spatiotemporal source analysis is 
the large spatiotemporal noise covariance matrix nec¬ 
essary for efficient estimators. This requires excessive 
amounts of data and can lead to unstable matrices [2], 
This makes it difficult to take temporal noise correla¬ 
tions into account. By Fourier transforming the sig¬ 
nals this problem can be largely overcome because the 
Fourier coefficients are assymptotically uncorrelated 


[3], so that only spatial noise covariances have to be 
estimated at each frequency. An additional advantage 
of analyzing the Fourier coefficients is that a subset of 
frequencies of interest can be selected, which may re¬ 
duce “noise” influences considerably. Frequency do¬ 
main dipole modelling has been developed in [5] and 
[ 6 ]. 

In this report we will assume a simple situation, in 
which there are two localized cortical areas that are 
well approximated by an equivalent current dipole 
source. Their electrical activities are related in that 
one of them is driving the other. Due to finite sig¬ 
nal transmission speed along neural axons, the signal 
from the driving source is leading in time with respect 
to the signal from the driven source. 


2 Methods 


Here we use a non-parameteric transfer function, in 
which a transfer coefficient for the transfer from one 
source to the other is estimated at each frequency: 


r,(k) = B(k)r,(k) rm 

-(SI’W'*, 


( 1 ) 


where A; is a frequency index and ff ]2 represents the 
transfer coefficient from the second source to the first 
at the fc-th frequency. The elements of rj represent 
the complex valued Fourier coefficients of the source 
amplitude signals, whereas the elements of £ repre¬ 
sent the innovations at the sources. We will refer to 
these equations as the structural model. 

The measurement model consists of the usual lead- 
field/amplitude specification, transformed into the 
frequency domain: 


Y(&) = C(0)rj(k) + e(k) (2) 
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Here, Y( k) is the complex valued vector of the sen¬ 
sor array measurements, 6 contains the dipole source 
locations and orientations, C is the leadfield matrix 



valued function, rj the source amplitudes as previ¬ 
ously defined, and e is the complex valued vector of 
the background noise activity. 

These equations imply the following cross-spectral 
matrix for Y(k): 

S(7 ,*) = 

C{6)[I-B{k)]r 1 C(ef + ® e (k) 

( 3 ) 

here |B(£;)}, |\I>(&)}, and {@ £ (A;)} are sequences 
of matrices for different frequencies k, and 7 = 
[0, (vec[B(&)]}, (vec[\I>(A;)]}, {vec [©(£;}}}]. The 
objective is to estimate 7 from the (multivariable) 
sample spectrum |S(A;)}. To this end the following 
negative loglikelihood function is minimized: 

L( 7 ) = ^{ln |S(7, k )| + tr[S(A;)E(7, *jf| 

k (4) 

- In \S(k)\ - (m + d)} 

where m is the number of sensors, and d the num¬ 
ber of dipoles. In general constraints have to be im¬ 
posed onB,^ and @ £ to make them identifiable. Of 
course there should be no more parameters than there 
are unique elements in (S(A;)}, but this is not a suffi¬ 
cient condition for identifiability. One class of models 
that is garanteed to be identified is the socalled recur¬ 
sive class, in which the matrix B can be forced into a 
lower triangular form by rearranging the order of the 
sources, and T' is diagonal [4], This is the situation 
that will be considered. 

3 Simulations 

A simulation was carried out to illustrate the method, 
and to assess the quallity of the estimates. 

In the simulation, the following time domain relation 
between two sources was introduced: 

>h(l) 0.m/ 2 (0'*Ki(0 

where (i(i) and (2 (t) are independent AR(2) pro¬ 
cesses (i.e. i &(k) is diagonal for all k). In the cur¬ 
rent simulation, the locations and orientations are as¬ 
sumed to be known, i.e. c(0) was fixed to its true 
value. In each of 500 simulations, a spectrum |S(A;)} 
was generated based on 32 timepoint epochs of con- 
tinous EEG. To compare the influence of sample size 
on the quality of the estimates, the generated spectra 
were based on 100 and 500 epochs of signal. 





Figure 1: estimated phases of the complex transfer 
coefficients at each frequency. The linear pattern is 
expected from the shift theorem in Fourier theory. The 
slope of the cross phase-spectrum is proportional to 
the time-lag that exists between signals. 


The estimates for the parameters in B, \I>. and @ £ 
were obtained by minimizing the log likelihood func¬ 
tion defined in equation 4. Table 1 summarizes the 
quality of the estimates. The parameters are approx¬ 
imately unbiased, and the quallity of the estimated 
standard errors is excelent. 

The well known shift-theorem in Fourier theory im¬ 
plies a linear relation between the phases of the trans¬ 
fer coefficients and the frequency index for the cur¬ 
rently simulated signals. Figure 1 shows this relation¬ 
ship where a box and whisker-plot of the estimated 
phases is given for 12 frequencies. The slope of this 
linear relation is proportional to the lag with which 
one signal is leading with respect to the other; the lag 
estimated from the phases averaged over simulations 
was -1.0015 time units, indicating that the amplitude 
signal of the second dipole {qf) was leading with re¬ 
spect to the first dipole (r/i ). In a second run, a model 
which assumed a transfer from the first dipole to the 
second (i.e. the other way around) was fitted to the 
simulated cross-spectral matrices. The resulting esti¬ 
mate of the lag was 0.9998 units, which is still indica¬ 
tive of an amplitude signal of the second source that 
is leading with respect to the first. 

4 Conclusion and Discussion 

Up till now interactions between different cortical ar¬ 
eas have only been studied by analyzing the raw sen¬ 
sor signals. Most common types of analysis that have 




Table 1: The quality of the parameter estimates. Bias: 
mean bias of the estimates of the nonfixed elements of 
either B(fe|, H?(k) and &(k). s.e.: mean estimated 
standard error of the parameters. s.eVs.d.: mean 
estimated standard error of the parameters devided 
the standard deviation of the the parameter estimates 
over 500 simulations respectively. 



100 epochs 

B © £ 

bias 

s.e. 

s.eVs.d. 

2.9e-4 2.0e-l 3.1e-3 

0.1 4.4e-l 9.8e-2 

1.0 9.9e-l 9.9e-l 


500 epochs 

B <l> @ £ 

bias 

s.e. 

s.eVs.d. 

5.1e-4 -5.7e-3 -1.8e-4 

0.03 0.2 0.03 

1.00 0.99 1.00 


been applied are coherence analysis and event related 
desynchronization. However these methods do not 
take into account the effects of volume conduction. 
These effects make it difficult to localize the activ¬ 
ity. Moreover they are also a major contributer to the 
correlations found between the sensor signals. There¬ 
fore we argue that the interactions between different 
brain areas are better analysed at the level of equiv¬ 
alent current dipoles than at the level of sensor sig¬ 
nals. The structural equations specification of the in¬ 
teraction between sources provides a framework for 
exploring and testing hypotheses concerning the cor¬ 
tical interactions. 

Although spatiotemporal source estimation is used 
widely in cognitive neuroscience research, the ob¬ 
tained estimates usually result from an ordinary least 
squares procedure. This procedure does not take into 
account the correlations between timepoints, and thus 
lead to inefficient estimators. Efficient estimators are 
obtained from a generalized least squares procedure, 
or maximum likelihood estimation. The reason that 
generalized least squares has not been applied is prob¬ 
ably because the estimation of the spatiotemporal co- 
variance matrix requires enormous amounts of data 
and can easily lead to unstable matrices. By trans¬ 
forming the EEG/MEG signal into the frequency do¬ 
main this problem is overcome in that the Fourier co¬ 
efficients are asymptotically uncorrelated. Therefore 
only the spatial correlation between spectral compo¬ 
nents need to be estimated and different frequencies 


can be analyzed separately, as is reflected in the sum¬ 
ming over k in eq. 4. 

The estimated quantaties in the frequency domain are 
usefull in their own right. The transfer function can 
be used directly to asses the amount of influence from 
one source to another. A lag estimate that might be 
obtained from the transfer coefficients can provide 
clues about the direction of the influence. 

One drawback of the Frequency domain analysis is 
that stationarity of the signals is assumed. When this 
assumption is not met by the data, the method may 
produce inaccurate results, and preferably other meth¬ 
ods have to be employed (e.g.[2j). 
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